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ABSTRACT 

We present a comparison of Spitzer IRS data for 51 OH megamaser (OHM) hosts and 15 non- 
masing ULIRGs. 10-25% of OHMs show evidence for the presence of an AGN, significantly lower 
than the estimated AGN fraction from previous optical and radio studies. Non-masing ULIRGs have 
a higher AGN fraction (50-95%) than OHMs, although some galaxies in both samples show evidence 
of co-existing starbursts and AGN. Radiative transfer models of the dust environment reveal that non- 
masing galaxies tend to have clumpy dust geometries commonly associated with AGN, while OHMs 
have deeper absorption consistent with a smooth, thick dust shell. Statistical analyses show that the 
major differences between masing and non-masing ULIRGs in the mid-IR relate to the optical depth 
and dust temperature, which we measure using the 9.7 fim silicate depth and 30-20 /im spectral slope 
from the IRS data. Dust temperatures of 40 — 80 K derived from the IRS data are consistent with 
predictions of OH pumping models and with a minimum Tdust required for maser production. The 
best-fit dust opacities (ry ~ 100 — 400), however, are nearly an order of magnitude larger than those 
predicted for OH inversion, and suggest that modifications to the model may be required. These 
diagnostics offer the first detailed test of an OHM pumping model based only on the properties of its 
host galaxy and provide important restrictions on the physical conditions relevant to OHM production. 
Subject headings: masers - galaxies: interactions - galaxies: nuclei - infrared: galaxies 



1. INTRODUCTION 

OH megamasers (OHMs) are 18-cm masers located in 
the nuclear regions of merging, (ultra)luminous infrared 
galaxies ([U]LIRGs). Possessing isotropic line luminosi- 
ties from 10 1 — 10 4 Lq, their hyperfine ratios, extremely 
broad linewidths, and large physical sizes point to a fun- 
damentally different o rigin t han the Galactic OH masers 
of the Milky Way (|Lol [2005). A rare phenomenon in the 
local universe (roughly 100 have been identified out to a 
redshift of z = 0.265), OHMs are exceptional probes of 
their environment due to their ability to be d etected at 
cosmic distances ([Darling fc Giovane lli 2002a). The as- 
sociation of the megamaser emission with merging galax- 
ies means that OHMs trace numerous extreme astrophys- 
ical processes, including high-intensity star formation, 
accretion in the central parts of galaxies, and the even- 
tual formation of massive black holes via binary black 
hole mergers. 

In order to use OHMs as tracers, however, the re- 
lationship between the maser emission and the envi- 
ronment of the host galaxies must be well quantified. 
Previous studies found no systematic difference between 
OH M hosts and ULIRGs of similar masses in the ra- 
dio (|Lonsdale et all 119981 : [Pihlstrom etafl I2005T ). opti- 
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cal (iBaan et al.l(l99l IDarling fc Giovanell [20061) . or X- 
ray ( Vignali et al.l I2005D regimes. OHM galaxies do, 
however, show exceptionally high dense gas fractions 
and h ave a distinctly non-linear IR-CO relation ([Darlingl 
12001 . Since OHMs are generated deep within the nu- 
clear regions of ULIRGs, however, the maser emission 
regions are almost always highly obscured, even at near- 
IR wavelengths. This means that observations capable of 
probing through the dust are critical both for determin- 
ing the parameters necessary for production of an OHM 
and determining its relation with the properties of the 
host galaxy. 

Mid-IR studies of OHM hosts to date are based pri- 
marily on photometry from the IRAS satellite; OHMs 
tend to occ ur in galaxies with color excesses at 25 
and 60 /zm (iHenkel et al.1 Il986[). high IR luminosities 
(|BaarJll989t IDarling fc Giovanell i 2002b), and steep far- 
IR spectral indices ([Chen et al.l I2007T ) . Spectroscopic 
studies of the mid-IR emission, however, offer much more 
powerful diagnostics that can explore the nature of the 
maser pumping mechanism and the associated OH emis- 
sion. We used the Infrared Spectrograph (IRS) aboard 
the Spitzer Space Telescope ([Werner et alJl20CM ) to ex- 
amine the nuclear regions of the merging OHM hosts. 
Mid-IR observations offer a particularly rich set of di- 
agnostics for ULIRGs, with measurements of AGN ac- 
tivity (high- ionization lines) , obscuring dust (absorption 
features and IR photometry), gas reservoirs (molecular 
absorption and H2 emission), and possible OH reservoirs 
(hydrocarbons and ices) all visible in the 5-35 fim region. 
For some galaxies, the masing gas can also be directly 
traced via the 34.6 /xm OH transition. 

We recently presented data from a spectroscopic sur- 
vey of 51 OHM hosts and 15 non-masing galaxies (Wil- 
lett et al. 2011; Paper I). In this follow-up analysis pa- 
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per, we present derived properties of the OHM galaxies 
and examine statistical differences between the two sam- 
ples. We also compare physical conditions in the mas- 
ing regions to those predicted from recent OHM pump- 
ing models. Finally, we describe how mid-IR diagnostics 
may serve as a useful selection technique for future OHM 
surveys. 

2. OBSERVATIONS AND DATA REDUCTION 

The Spitzer data for both the OHMs and the non- 
masing galaxies come from multiple observing programs; 
we observed 24 OHMs in a dedicated Cycle 3 program 
(30407) for IRS observations of OHM hosts. Additional 
data for both OHMs and confirmed non-masing galaxies 
were drawn from the Spitzer archive, with approximately 
half from the IRS GTO sample of ULIRGs. In order to 
create a uniform data set, we required that all galaxies 
have full coverage in both the low- and high-resolution 
modules. Including objects from our dedicated OHM 
observing program, the samples contain 51 OHM hosts 
and 15 non-masing galaxies with no OH emission above 
Loh = 10 23 L Q . 

While OH observations cover the majority of radio 
galaxies and ULIRGs in the local universe (z < 0.2) 
which are likely candidates for OHM emission, our sam- 
ple was largely constrained by the availability of data 
from the Spitzer archive. This sample is not complete, 
as there are many objects (both OHMs and non-masing 
galaxies with firm upper limits) for which mid-IR spec- 
troscopy was not available. 

Since the archival objects did not come from a uni- 
fied observing program, the version of the Spitzer data 
pipeline and the level of processing vary slightly from ob- 
ject to object - we used the most recent versions available 
in the archive (vl5.3.0 or later). The reduction pipeline 
is described in detail in Paper I; briefly, we use the basic 
calibrated data (BCD) products, subtracting background 
sky for all LR modules and medianing subsequent ex- 
posures to remove transient effects. The images were 
cleaned using the IDL routine IRS CLE AN_M ASK and 
the 1-D spectra then extracted using the Spitzer IRS 
Custom Extractor (SPICE) v. 2.0. Almost all galaxies 
in our sample were unresolved with the IRS and were 
treated as point sources. 

Low-resolution modules were stitched together to 
match continuum levels by using a multiplicative scaling, 
fixing the LL1 module and then stitching the other three 
modules at the points of overlap. The spectra were cali- 
brated as a single unit to 22 or 25 ^m photometry. Noisy 
regions (typically 10-30 pixels) corresponding to areas of 
lowered detector sensitivity at the edges of the SH and 
LH orders were trimmed from the final 1-D spectra. In 
isolated cases, exceptional 1-channel features appearing 
in only a single nod were either manually removed or re- 
placed using the data from the uncorrupted nod position. 

3. DISCUSSION AND ANALYSIS 

3.1. Spectral energy distributions 

The low-resolution (LR) IRS spectra are powerful in- 
dicators of the overall spectral shape of the galaxies, 
typically dominated by reprocessed emission from dust 
heated by star formation and /or an active nucleus (e.g., 
lArmus et aLll2007t lHao et al.ll2007|) . Figure []] shows the 
median LR spectra for all galaxies in both the OHM and 
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Figure 1. Medianed low-resolution spectra for all OHMs (black) 
and non-masing galaxies (red). The lcr error bars for each medi- 
aned pixel are also shown. The dashed line shows where the fluxes 
are normalized at 15 /im. 

non-masing samples. A clear difference in the spectral 
shape between the two samples is apparent; the OHMs 
show deeper absorption at both 9.7 and 18 fim and 
steeper continuum from 15-35 /im. Discrepancies in in- 
dividual emission and absorption features are also ap- 
parent; the PAH emission at 7.7 /im is broader in the 
OHM template, with the 8.6 /im feature largely sup- 
pressed (possibly due to extinction from silicate dust). 
Similarly, the H2 S(3) A9.67 line is clearly seen in the me- 
dian OHM template and suppressed in the non-masing 
sample. 

The medianed OHM spectra also reveals a clear ab- 
sorption feature n e ar 6 associated with water ice 
HSpoon et alJ[2002L l200l : the same feature is not seen 
in the medianed template of the non-masing galaxies. 
This is consistent with individual detection rates in the 
two samples (24/51 OHMs, 3/15 non-masing galaxies). 
Since water ice is a possible reservoir for the masing 
OH molecules in their ISM gas-phase, distinct differ- 
ence between the two populations have implications for 
OHM emission. If large fractions of the available OH 
are locked up in solid forms (ice mantles on dust grains, 
for example), then the reservoir of gas-phase OH could 
be depleted to a degree that would quench maser emis- 
sion. This could be due to a harder radiation environ- 
ment in the non-masing ULIRGs, as sufficiently strong 
UV radiation can dissociate OH even in the ice phase 
(|Andersson fc van Dishoeckl [20081) . 

3.2. Narrow-line region gas 

We traced the hardness of the radiation field by com- 
paring the excitation states of the fine-structure neon and 
sulfur lines, plotting the ratio [Ne III]/[Ne II] against 
the [S IV]A10.5/[S III] A18.71 (Figure©. Detection of 
all four lines occurred in less than 50% of the sample 
(14/51 OHMs, 8/15 non-masing galaxies); non-detection 
of [S IV] is the limiting factor for almost all galaxies. 
Since [S IV] lies near the silicate absorption at 9.7 /xm, 
extinction caused by mixing of the dust and ionized gas 
may suppress observation of this line for dust-rich galax- 
ies. 

The line ratios from our samples are compared to 
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Figure 2. Diagnostic of the excitation state in the narrow-line 
region for OHMs (red) and non-masing galaxies (blue) in which 
all four of the [Ne II], [Ne III], [S III], and [S IV] t ransitions are 
detec ted. Also shown are line ratios fo r ULIRGs l|Farra h et alj 
[2007D , active galaxies ISturm et al.lp"0 02; Tomr najm et aI.II200aT T 
and starburst galaxies (Vcrma et al. 2003). 

larger popula tions of ULIRGs (iFarrah et al.ll2007f). ac - 
tive galaxies (ISturm et al.ll2002t iTommasin et al.ll2008f l. 
and starbursts (|Verma et al.ll20031 ). All galaxies show 
a correlation between higher ionization states for both 
species, with galaxy types relatively evenly distributed 
through the total range of line ratios. Fits for each set 
are consistent with a slope between 0.5 and 1.0; although 
OHMs have a shallower slope (0.5 ± 0.7) than the non- 
masing galaxies (0.8 ± 0.7), the uncertainties in both fits 
are too large to distinguish them from each other or the 
larger population of ULIR Gs. These slopes are also con- 
sistent with the results of IDale et al.1 (|2006D . who found 
that nuclear regions of galaxies in the SINGS sample ex- 
hibited a similar trend (although they used [S III] A33.48 
instead of [S III] A18.71). Mixing of the line ratios for 
OHM hosts and non-masing galaxies, along with the lack 
of a clear locus for either population, suggests that the 
ionization state of the narrow-line gas is not a factor in 
triggering an OHM. 

4. DERIVED PROPERTIES 
4.1. Velocities 

The IRS high-resolution spectra contain multiple nar- 
row lines that can be used for accurate redshift mea- 
surements. We computed a systemic IR velocity from 
the weighted mean of all detected HR line centroids, 
typically fixed by the strongest transitions (including 
H 2 , [Ne II], and [Ne III]). The scatter of the individual 
species around the mean velocity is ~ 200 km/s for the 
strongest lines, while the velocity resolution of the IRS 
is ~ 500 km/s. There is also a possible unknown veloc- 
ity component from the Spitzer spacecraft, which may 
be as high as 30 km/s. We found no statistically signif- 
icant trend for individual species versus redshift, similar 
to the trend found in OHM host optical line redshifts by 
iDarling fe Giovanellil (|2006f ). 

We compared both the individual and systemic IR ve- 
locities to those measured from opt ical spectroscopy and 
to the velocity of the OHM itself. IDarling fc Giovanellil 
( 2006) found a significant asymmetry in the OHM-optical 
redshift distribution, with the OHM emission somewhat 



Figure 3. OH-IR velocity offset vs. the systemic IR velocity for 
the OHMs. The mean velocity offset (dashed) for the sample is 
—90 ± 19 km/s, showing a blueshift of the OHM relative to the 
IR emission lines. Plotted uncertainties are statistical only, and do 
not account for the unknown velocity (< 30 km/s) of the spacecraft 
at the time of observation. 

blueshifted with respect to the optical emission. The re- 
sults for the mean OH-IR velocity offset show a similar 
blueshift of Av avg = -90 ± 19 km/s (Figure EJ). This 
is consistent with our measurement of no systematic off- 
set between the IR and optical velocities {Av op t-iR — 
— 13 ± 20 km/ s). While this agrees with the results of 
iDarlind (|2007D . the optical/IR agreement is somewhat 
puzzling given the large amounts of dust in the actively 
merging galaxies. If the IR lines truly come from the 
nuclear regions and the optical lines from superficial gas, 
an offset between the two sets of transitions might be ex- 
pected - non-detection of this effect (and relatively small 
scatter) may imply that many of the IR lines are super- 
ficial. This is supported by the detection of H 2 S(3) and 
[S IV] emission on top of the 9.7 /xm dust absorption fea- 
ture. We note that for high-ionization lines that must 
originate near the nucleus ([Ne V] 14 and 24 /mi), we 
do not have enough detections to measure a significant 
statistical offset. 

The alignment of the mean IR and optical veloc- 
ities could be partly due to a selection effect, since 
the lines arc primarily identified on the basis of pre- 
existing optical redshifts (although mis-identification of 
lines wo uld require off s ets o f thousands of km/s or 
greater). iSpoon fc Holtl (|2009l ) showed that [Ne III] and 
[Ne V] emission in ULIRGs can be offset by more than 
200 km/s, likely explained by decelerating outflows that 
are photoionized by AGN. The lack of a systematic 
blueshift of the OHM in our sample may be an indi- 
cator that outflows are not common in the host galaxies, 
a further indication of a tendency for OHMs not to be 
associated with AGN. A two-sided Kolmogorov-Smirnov 
test showed no significant difference in Av op t-iR for the 
O HM and non-masing p opul ations. 

IDarling fc Giovanellil (|2006| ) also found a weak corre- 
lation between the magnitude of the OHM blueshift and 
the strength of the OHM (as measured by log Lqh and 
the linewidth W W67 ). The blueshift of the OHM with re- 
spect to the IR emission showed no significant correlation 
for either parameter for the galaxies in our sample. 
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Figure 4. Luminosity of [Nc II] + [Nc III] lines as a function 
of Lin for t ne OHM (diamond) and non-masing (cross) samples. 
Linear fits to both samples are shown by the solid (OHMs) and 
dashed (non-masing) lines; the d o tted line shows the fit for the 
much larger sample of Ho k Ketol 020071 ). The fit to the OHMs is 
within the scatter of both the other two samples; fits to the HK07 
and non-masing galaxies, however, do not agree. 

4.2. Star formation 

We also examined the relationship between the OHM 
and the star formatio n rate (SFR) in the host galax- 
ies. iHo fc Ketol (|2007t l use the fine-structure [Ne II] and 
[Ne III] lines as diagnostics in galaxies spanning more 
than five decades of IR luminosity. Neon emission is 
a useful tracer for the SFR due to its abundance in 
HII regions, ionization energies that make the singly- and 
doubly-ionized species among the primary coolants for 
gas heated by massive stars, and relative insensitivity to 
dust extinction (particularly when compared to common 
tracers in the optical/UV such as [O III] A5007A and 
Ha). 

Figure [4] shows the relationship between the far- 
infr ared luminosity (Ljr\ m easured using the method 
of iSanders fc Mirabel Il996f ) and the integrated lumi- 
nosity of the neon lines. We plot the results for both 
samples in Figure 0] along with a least-squares linear fit. 
We also o verplot the relatio ns found from the broader 
sample in lHo fc Ketol ((2007). Although there is a mod- 
erate correlation between L^ e and Lm (as is expected 
for any comparison involving two luminosities) , the scat- 
ter is considerable. The combined neon luminosities for 
the OHMs yield a fit of 



log [L NeII+NeIII ] = (1.0 ± 0.5) log [L IR ] - (4 ± 6), (1) 

with both luminosities measured in L@. The fit to the 
non-masing galaxies is: 



log [L NeII+NeIII ] = (1.5 ± 0.4) log [Lin] - (10 ± 5). (2) 

Both the slope and offset for the OHMs are consis- 
tent with the relationship found by IHo fe Ketol (|2007j ) : 
log [L NeII+NeIII ] = (0.98 ± 0.069) log [L IR ] - (2.78 ± 
0.70). The fits for the OHMs and non-masing galaxies 
are also consistent w ithin the large scatter. The fits to 
the lHo fe Ketol (l200l and non-masing galaxies, however, 
do not overlap (within their respective lcr scatter), indi- 



cating a possible marginal difference in SFR. 

The larger uncertainties in the OHM and non-masing 
galaxies' slopes are attributed to their narrow range in 
Liu. Both samples have a lower limit on Lm that lies at 
the high end of the IHo fe Ketol (|2007n data. The upper 
end of the Lm range reflects the low space density of 
HyLIRGs with L m > 10 12 ' 5 L Q . The net effect yields 
only ~ 1.5 dex of Lm over wh i ch a r elation can be fit; 
since the galaxies in IHo fe Ketol (2007) is over more than 
five decades of Lm , their correlation is much tighter. 
iFarrah et a l. (2007) suggest the offset between the slopes 
is a result of higher extinction in the nuclear regions of 
ULIRGs relative to lower-luminosity starbursts that fix 
the height of the IHo fc Ketol (|2007f) relation. The fact 
that both samples are neon-underluminous compared 
to the larger data set agrees with the high extinction 
(ry ~ 300) found by fitting dust models ( §4.50 . Since the 
non-masing galaxies have even lower neon fluxes, their to- 
tal Lm likely has a lower overall contribution from star 
formation. 

We computed star formatio n rates for our sam ples us- 
ing two diagnostics: L^ e from IHo fc Ketol (l2007f) and the 
starbu rst far-IR luminosity calibration from iKennicuttl 
(1998) (Table Q}. For the neon relation, we assume an 
ionization fraction of fi on — 0.6 and neon ionization frac- 
tions of fNe+ = 0-75 and fNe++ = 0.1- The OHM mean 
star formation rate is (SFR]y e ) = 120 ± 16 M /yr, com- 
pared to (SFR.N e ) = 65 ± 21 M /yr for non-masing 
galaxies. Using the Lpm, calibration with IRAS pho- 
tometry yields a higher average SFR and larger scatter 
for both samples, with (SFR FIR ) = 300 ±30 M Q /yr for 
OHMs and {SFR Ne } = 210 ± 50 M /yr for non-masing 
galaxies. The closer agreement between the two popula- 
tions for this calibrator reflects the Lm criterion used in 
selecting the non-masing sample. 

Both the large scatter and the difference in SFR be- 
tween the diagnostics illustrate the difficulties in charac- 
terizing a local phenomenon over a large volume. This 
emphasizes the fact that all OHMs are mergers with mul- 
tiple components, with sites of star formation likely sep- 
arated by tens of kpc. The linear correlation between the 
two SFRs is also relatively low, with a Spearman's rho 
of p = 0.54. This could indicate a component for heating 
the dust that does not come from star formation, such 
as an AGN. 

4.3. AGN vs. starburst 

A key issue surrounding ULIRGs is their central source 
of power - does it come from AGN or starbursts? 
iBaan et ahl (|1998f ) use optical classifications of OHM 
hosts to claim that 45% of the host galaxies show signs 
of a pure AGN (Seyfert and LINER spectra), with 
an additional 22.5% displaying composite spectra with 
characteristics from both A GN and starburst activity. 
iDarling fc Giovanellil (2006) compare a sample of OHM 
host galaxies vs. non-masing ULIRGs and find that 
42% are LINERs, 25% Seyfert 2 galaxies, and 33% star- 
bursts; classifications are similar for both samples. They 
also find few significant correlations between the OHM 
emission and the optical properties of their host galax- 
ies. Classification using the radio and FIR properties of 
the nuclei, however, show only 34% of the sa mple with 
AGN characteristics (|Baan fc Klockner] 12006): these in- 
clude multiple objects optically classified as LINERs or 
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Table 1 

Derived mid-IR and radio properties for OHMs and non-masing galaxies 
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1.5 


56 


41 


-0.6 


-3.0 


0.0 


IRAS 13428+5608 


110 


244 


65 


250 


0.0 


200 


67 


-0.4 


-1.5 


-1.5 


IRAS 13451+1232 


174 


253 


68 


1000 


0.0 


31 


14 


-0.0003 


-2.0 


0.0 


IRAS 14059+2000 


104 


149 


63 


950 


0.0 


61 


43 


-1.1 


-0.5 


0.0 


IRAS 14070+0525 


341 


1092 


95 


250 


0.0 


410 


64 


-1.1 


-1.5 


** 


IRAS 14553+1245 


122 


128 


73 


250 


0.0 


190 


68 


-0.6 


-2.0 


-1.5 


IRAS 15327+2340 


24 


266 


60 


250 


0.0 


410 


64 


-0.6 


-1.5 


** 


IRAS 16090-0139 


200 


607 


62 


300 


0.0 


190 


63 


-0.4 


-1.5 


-1.5 


IRAS 16255+2801 


64 


151 


70 


250 


0.0 


210 


67 


-0.9 


-2.0 


-1.5 


IRAS 16300+1558 


223 


927 


54 


300 


0.0 


410 


59 


-0.3 


-1.5 


** 


IRAS 17207-0014 


94 


456 


59 


300 


0.0 


410 


59 


-1.0 


-1.5 


** 


IRAS 18368+3549 


129 


299 


56 


300 


0.0 


260 


61 


-0.2 


0.0 


-1.0 


IRAS 18588+3517 


94 


144 


71 


250 


0.0 


280 


66 


-0.8 


-2.0 


-1.0 


IRAS 20100-4156 


173 


732 


62 


250 


0.0 


410 


64 




-1.0 


** 


IRAS 20286+1846 


16 


201 


37 


250 


0.0 


230 


67 


-1.4 


0.0 


-1.5 


IRAS 21077+3358 


167 


177 


76 


200 


0.0 


410 


70 


-0.4 


-2.5 


** 


IRAS 21272+2514 


75 


151 


77 


200 


0.0 


410 


70 


-1.5 


-2.0 


** 


IRAS 22055+3024 


109 


267 


71 


200 


0.0 


240 


73 


-0.7 


-2.0 


-1.5 


IRAS 22116+0437 


178 


225 


76 


200 


0.0 


380 


71 


-0.2 


-2.5 


** 


IRAS 22491-1808 


44 


246 


58 


250 


0.0 


410 


64 


-1.1 


-1.5 


** 


IRAS 23028+0725 


79 


125 


91 
















IRAS 23233+0946 


118 


232 


67 


200 


0.0 


380 


71 


-0.3 


0.0 


** 


IRAS 23365+3604 


43 


245 


70 


200 


0.0 


410 


70 




-2.0 


** 


Non-masing IRAS 00163—1039 


69 


40 


52 


1000 


0.0 


50 


43 




0.0 


0.0 


IRAS 01572+0009 


559 


513 


81 


850 


0.0 


23 


■18 




—0.5 


0.0 


IRAS 05083+7936 


181 


148 


43 


1000 


0.0 


34 


44 




0.0 


0.0 


IRAS 06538+4628 


26 


30 


51 


75 


0.0 


380 


109 




0.0 


** 


IRAS 08559+1053 


288 


261 


74 


1000 


0.5 


46 


42 




-1.5 


0.0 


IRAS 09437+0317 


4 


24 


38 


1000 


0.0 


16 


43 




0.0 


0.0 


IRAS 10565+2448 


139 


188 


55 


950 


0.0 


67 


43 




-0.5 


0.0 


IRAS 11119+3257 


196 


519 


78 


150 


1.0 


51 


82 




-2.0 


-2.5 


IRAS 13349+2438 


69 


42 


21 1 


950 


2.0 


90 


40 






0.0 


IRAS 15001+1433 


318 


454 


67 


1000 


0.0 


55 


43 




-1.5 


0.0 


IRAS 15206+3342 


591 


237 


78 


1000 


0.0 


31 


44 




-1.5 


0.0 


IRAS 20460+1925 




186 


112 
















IRAS 23007+0836 


50 


45 


79 


650 


0.0 


28 


53 




-1.5 


0.0 


IRAS 23394-0353 


28 


22 


11 


700 


0.0 


56 


49 




0.0 


0.0 


IRAS 23498+2423 


662 


477 


93 


500 


1.0 


45 


53 




-3.0 


0.0 



Note. — ^du^t ls dust temperature fit to the photometric greybody; the second Tdust is the temperature at the outer envelope 
fit by the DUSTY models. t"^P 7 is the apparent mascr optical depth calculated from radio fluxes in the literature (Equation^, while 
r i667 anc ^ T i66'f TY are predicted values from the LE08 model based on the data in £|6. 11 and j|6.21 respectively. ** indicates that the 
predicted T1667 fell outside the available contours for the LE08 model. 
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OHMs 


[Nc V] 
— no [Ne V] 






Non-masing 



5678910 15 20 30 

Wovelength [/u.m] 
Figure 5. IRS low-resolution spectra, showing the difference 
in the median-stacked spectra for galaxies with [Ne V] emission 
(black) and with [Ne V] upper limits (red). [Ne V] was detected in 
4/51 OHMs (top) and 8/15 non-masing galaxies (bottom). Spectra 
are normalized in flux at A = 15 fim (dotted line). 

composite objects that show no AGN activity in the ra- 
dio. It is suggested that the differences in classification 
lie in the large amounts of extinction at optical wave- 
lengths due to dust obscuring the nucleus. 

In the mid-IR, high-ionization fine structure emission 
lines are the simplest and most unambiguous tracers of 
AGN activity. [Ne V] has an ionization energy of 97.1 eV, 
a level typically too high to be reached by young O 
and B stars. [O IV] has a smaller ionization energy 
of 55 eV, which is often seen in AGN and in several 
optically-identified starburst galaxies. In contrast, the 
[Ne V] line is only seen in integrated galactic spectra that 
harbo r AGN, although it is not ubiquitous - iDudik et al.l 
(|2007h detect [Ne V] in 19/41 Seyfert and LINER galax- 
ies, for example. It is possible in these cases that lines 
are present, but that differential extinction in the mid- 
IR obscures their emission. Since AGN occupy a much 
smaller volume (< 1 pc) than a typical starburst and 
have harder radiation fields, the high-ionization regions 
where the neon and oxygen are emitted are likely to be 
more deeply obscured than their low-ionization counter- 
parts. 

[Ne V] A14.3 is detected in 4/51 OHMs and 8/15 non- 
m asing galaxies (Pap er I) . ULIRGs in the larger sample 
of lFarrah etlll (|2007t ) show [Ne V] A14.3 in 22/53 galax- 
ies, three of which are OHMs that overlap with our sam- 
ple. [O IV] is seen in 21/22 galaxies in their sample that 
display [Ne V] and only in two that do not, demonstrat- 
ing a close but not perfect association. The difference 
in detection rates suggests that the presence of an OHM 
selects against AGN with high-ionization emission; this 
may relate to the timescale of the galactic merger. If 
OHMs are associated with a particular phase since the 
onset of the host galaxies' merger (and possibly a de- 
lay before the activation of the AGN), then this would 
explain why the OHM sample has so few [Ne V] detec- 
tions compared to non-masing galaxies and ULIRGs in 
general. 

Given the strong differences in the [Ne V] detection 
rate between the samples, we examined whether galax- 
ies emitting [Ne V] might reveal other parameters rele- 



vant to OHM formation. In the average IRS LR spec- 
tra for OHMs (Fig. [5]), galaxies with [Ne V] emission 
show a shallower 30-20 fim slope than galaxies without 
high-ionization lines (030-20 = 3.7 vs. 5.4). The sili- 
cate depths and PAH luminosities, however, are broadly 
consistent for both samples. For non-masing galaxies, 
Figure [5] shows that the average 9.7 /xm silicate depth is 
shallower for galaxies that show [Ne V] (6*9.7 = 0.5 vs 
0.7), but that the 030-20 slopes are similar. The EW 
of the 6.2 /im PAH feature is also smaller in galaxies 
with [Ne V] ; the high-ionization lines are consistent with 
the presence of an AGN that dissociates large molecules. 
For the OHMs, the presence of [Ne V] shows no effect on 
either Loh or the peak flux at 1667 MHz. 

Other mid-IR diagnostics can also be used to charac- 
terize the contribut ion of AGN and/or starburst features. 
ISpoon et afl (j2007t ) plot the PAH 6.2 fim EW against 
the silicate 9.7 /j,m strength in a "fork" diagram. IRS 
data show two distinct branches of galaxies for this diag- 
nostic, with one representing a largely AGN-dominated 
population (weak PAH emission and little to no silicate 
absorption) and another containing ULIRGs/HyLIRGs, 
obscured AGN, and starburst galaxies (stronger PAH 
emission coupled with deeper silicate stren gths). 

We reproduce the fork diagram from ISpoon et aLl 
(2007) with our IRS data overlaid on the broader sam- 
ple of ULIRGs, starburst galaxies, and AGN in Figure [6j 
OHMs lie almost exclusively along the top branch and 
share significant overlap with optically identified star- 
burst galaxies, which typically have strong PAH emission 
but weak to moderate silicate absorption. The locus of 
the OHMs on the fork diagram agrees with the [Ne V] and 
[O IV] data; only four OHMs lie on the horizontal, AGN- 
dominated branch. The non-masing galaxies are prin- 
cipally found along the horizontal branch, with a wide 
range of PAH EW but lower S 9 . 7 than the OHMs. A 
small region of overlap does exist between the two sam- 
ples near the "knee" (high PAH EW and weaker silicate 
absorption). The absence of non-masing ULIRGs on the 
upper branch is one of the first clear spectral diagnos- 
tics of OHMs, based only on the properties of the host 
galaxy. 

Using the mid-infrared diagnostics, we estimate the 
AGN contribution by assuming that [Ne V] clearly in- 
dicates an AGN and that [O IV] detection or placement 
on the horizontal branch in the fork diagram indicates a 
possible AGN. Based on the IRS data, the AGN fraction 
of OHMs is between 10 and 25%, compared to a much 
higher fraction of 50-95% for non-masing ULIRGs. The 
AGN contributions for the combined samples are con- 
sistent with that estima ted from radio/FIR diagnostics 
(jBaan fc Klockne"rll2006[ ). The co-existence of AGN and 
starbursts in some nuclei is also supported by the pres- 
ence of galaxies showing both high-ionization emission 
and large PAH equivalent widths. The mid-IR AGN frac- 
tion for OHMs, however, is significantly low er than that 
estim ated from optical diagnostics (45-70%; iBaan et al.l 
[1991 . which can be significantly affected by dust obscu- 
ration around the nuclear regions. 

4.4. Dust temperatures 

We computed dust temperatures for both samples us- 
ing two methods. The first is a single fit to the in- 
tegrated IR emission by assuming a single-temperature 
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Figure 6. "Fork" diagram displaying the 6.2 /im PAH equivalent 
width vs. the silicate strength at 9.7 fim. OHM galaxies from 
both our program and the archive are shown in red; non-masing 
galaxies are shown in blue. Top: binned distribution of PAH EW 
for both samples; Ri ght: binned d i stribu tion of S9.7. Additional 
Spitzer data are from lSpoon et alj ll2007t) . 

modified blackbody as a template (the second is taken 
from the model fits i n §4.5p . We adapt the broad SED of 
lYun fc Carilll (l200l . where the emission follows a ther- 
mal blackbody above a critical frequency v c where the 
dust clouds become optically thick and a greybody spec- 
trum below v c . For an object subtending an angular di- 
ameter [arcsec], the expected flux density at frequency 
v [GHz] is 



S d [v] = 2.8 x 10" 



e 0.048iVT d - 1 



1 



Jy- (3) 



In our analysis, 6 is a free parameter (accounting 
for the dependence of flux density on distance) and as- 
sumed v c — 2000 GHz (150 fim) and an emissivity index 
ft = 1.35. The fits are relatively insensitive to the choice 
of P since for most galaxies we lack photometric data 
points below the critical frequency. We used fluxes from 
IRAS and IRS peakups to fit the curve with photometry 
from 12-100 /im, allowing both the physical temperature 
and the peak intensity (a function of both distance and 
extinction) to vary. The majority of the galaxies only 
have IRAS detections at 60 and/or 100 /im, in addition 
to the IRS peakups. 

The mean temperature for the OHMs is Tdust = 
66 ± 12 K, while the non-masing galaxies have T^ ust = 
80 ± 50 K. These uncertainties are the statistical la en- 
velopes for the sample and do not address the physical 
relevance of fitting the galaxies with a single tempera- 
ture fit. The hottest temperature measured is in the non- 
masing galaxy IRAS 13349+2438 (T dust = 243 K), which 
is more than twice as hot as the next-highest galaxy. The 
greybody temperatures for all galaxies are listed as T% U8t 
in Table [U 

4.5. Modeling the dust environment 

The infrared emission of ULIRGs is dominated by ra- 
diation from heated dust; thus, differences in the distri- 
bution of dust have significant influences on the mid-IR 
spectra. Since the mid-IR photons are responsible for 




Figure 7. Feature-feature diagram plotting the relative strengths 
of the 9.7 and 18 /im silicate features. The lines represent 
models of different d ust geometries for cool, oxygen-rich silicates 
HSirockv et al. 2008). Dotted, dashed, and solid lines represent 
the thickness of the dust shell in the smooth models (Y = 
Router I Rinner) with a flat radial density distribution (q = 0, 
where p[r] oc r~ q ). The green tracks model clumpy geometries 
with varying numbers of dust clouds located along the line of sight 
(TVrj = 1,3, 5 from upper right to lower left). The black star in the 
upper right corner is the starting point for all optically thin dust in 
the models (59.7 = 1.26, Sis = 0.67). Silicate strengths from the 
IRS data are shown for the OHMS (red) and non-masing galaxies 
(blue). 

maser pumping, this may also have an important effect 
on the presence of OHMs in ULIRGs. Using the radia- 
tive transfer code DUSTY, we modeled the dust environ- 
ment of the galaxies in our sample for two geometries: a 
smooth, thick dust shell and a clumpy torus. 

The DUSTY code models the dust environment as 
a smooth, s pherical distribution of centrally illumi- 
nated dust (jlvezic fc Elitzurl I1997D . Our models as- 
sumed the_dust^_wascom£osed of cool, oxygen-rich sili- 
cates (Ossenko pf et al.lll992D and that the heating source 
follows a broken-power law luminosity function. The 
code fits for the thickness of the dust shell (Y = 
Router I Rinner) , the power-law index q of the radial den- 
sity profile (p[r] oc r~ q ), and the total optical depth Ty 
at 0.55 fim. The inner dust radius, which is defined by 
the source luminosity and dust sublimation temperature, 
is the only free physical parameter in the code. DUSTY 
then generates a grid of artificial spectra at a variety of 
radii, from which parameters such as the dust tempera- 
ture can be extracted. 

Motivated by the evidence that some fraction of our 
galaxies host AGN, we generated a second set of mod- 
els for a clumpy distribution of dust, which may bet- 
ter represent the environment around active galaxies 
(jLandt et all 120101 ) . We calculate the source function 
for individ ual clumps using DUSTY a nd used the code 
CLUMPY (INenkova et alll2002L l2008aTbh to account for 
the new geometry. CLUMPY assumes a distribution of 
individual dusty clouds in a torus around the central il- 
luminating source. The code fits for Y, q, and Ty as 
well as the number of clouds along the line of sight (No), 
the angular dependence of cloud distribution away from 
the equatorial plane (a), and the inclination angle of the 
galaxy (i). Neither a nor i were well-constrained param- 
eters in our models. 
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We first used both the DUSTY and CLUMPY models 
to examine the overall dust distributions of the galaxies 
in our sample. The "featu re-feature" diagram, devel- 
oped bv lSirockv et al.1 (12008(1 . plots the depths of the 9.7 
and 18 /im silicate features against each other. These 
are compared to tracks from the radiative transfer mod- 
els; we used a small set of dust geometries and plot- 
ted the expected silicate ratios for a large ra nge of op- 
tical depths. Following iSirockv et al.l (J2008), we gen- 
erated tracks for three smooth geometries that vary in 
shell thickness (Y = 100, 200, 400 for q = 0.0) and three 
clumpy geometries that vary in the number of clouds 
(jVo = 1,3, 5 for q = 0, Y = 30); the optical depth is 
then allowed to vary for each model from to 80. 

Figure [7] shows the tracks for the different dust ge- 
ometries, as well as the measured silicate ratios for the 
OHMs and non-masing galaxies from the IRS spectra. 
The non-masing galaxies occupy a much smaller locus 
of possible dust geometries than the OHMs, showing no 
deep absorption (S < —1.2) in either silicate feature. As 
a result, most non-masing galaxies are best fit by one of 
the clumpy dust geometries. The OHMs occupy a much 
larger region; while a few galaxies fall close to the clumpy 
tracks, the majority of OHMs have deep 9.7 /xm absorp- 
tion that only be achieved wit h a smooth, embedding 
medium. iLevenson et al.l (|2007l) show that such deep ab- 
sorption requires a large temperature gradient across the 
absorbing medium, which can only be achieved if the dust 
screen is both geometrically and optically thick. While 
the silicate ratios are not sensitive enough to strongly 
constrain either Y or No, it does demonstrate a clear 
difference in the dust environments of the two ULIRG 
populations. 

Puzzlingly, most of the OHMs fall below the tracks pre- 
dicted for the smooth dust geometries at S9.7 < —1.5. 
If this is a systematic effect, then this implies that ei- 
ther the 18 /im feature is being overestimated or the 
9.7 /im feature underestimated with respect to the mod- 
els. We consider the former more physically probable; if 
the OHMs were shifted to the left in Figure[7]to lie on the 
smooth tracks, this would imply absorption depths of up 
to S9.7 ~ 6 — 7, much deeper than any seen in a ULIRG 
to date. As a check, we directly co mpared our measured 
silicate depths to those published in ISirockv et al.l (2008) 
for ten galaxies that appear in both samples; however, no 
systematic difference in absorption strengths was found. 
Assuming the methods are consistent, this may indicate 
that a different geometry must be implemented in the 
code for the most heavily embedded ULIRGs. 

Figure [7] suggests that OHMs are best modeled by a 
smooth geometry; based on these results, we attempted 
to further constrain the dust geometry for our galaxies by 
fitting each IRS spectra with the smooth-shell model. We 
used DUSTY to generate a grid of 13860 artificial spec- 
tra; the parameter ranges are in Table O We chose these 
values to span the expected physical range for ULIRGs: 
these include shell thicknesses out to Y = 1000 (400 pc 
for typical values of dust sublimation temperature and 
the heating source luminosity), power-law indices from 
— 3, and Ty extending up to 500. 

After generating the grid of artificial spectra, we 
needed to identify the best fit for each IRS spectrum. 
Since DUSTY only models continuum and dust features, 
we removed the PAH, atomic, and molecular lines from 



Table 2 

Grid parameters for DUSTY model fits 



DUSTY grid 


Y 


2, 5, 10, 15, 20, 30, 50, 75, 100, 150, 200, 300, 




400, 500, 750, 1000 


g 


0, 0.5, 1.0, 1.5, 2.0, 2.5, 3.0 




0.1 - 500 



Note. — rv is binned on a logarithmic scale with 180 
steps between 0.1 and 500. 



the IRS data to improve the quality of the fit. For this 
we employed PAHFIT, a set of IDL routines that per- 
forms spectral deco mposition of low-resolution IRS data 
(ISmith et all 12001 . While removing narrow line emis- 
sion was typically clean, subtraction of the PAH emission 
often increased the area of the 9.7 /j,m silicate feature, 
since the wings of the 8.6 and 11.3 PAH profiles fill in the 
dust absorption. Once the IRS spectra were reduced to 
continuum + dust features, we re-binned the data to the 
resolution of the DUSTY grid and found t he best fit fol- 
lowing the error minimization technique of Nikutta et al.l 
( 2009). Results of the best fit Y, q, Ty, and Tdust at the 
outer edge for each galaxy are given in Table [TJ 

Figure [5] shows the distributions of the model param- 
eters for the best DUSTY fits to the OHMs and non- 
masing galaxies. The best fits have a uniformly flat den- 
sity profile for almost all galaxies in our sample; only 
3/51 OHMs and 4/15 non-masing galaxies had best fits 
with q > 0. The two samples also have similar dust tem- 
peratures, with OHMs slightly warmer on average than 
the non-masing ga laxies ((Tdust) = 62 vs . 53 K), confirm- 
ing the results of iDarling fc Giovanellil (|2002al ). These 
values are consistent with the greybody dust tempera- 
ture measured with IR photometry (Equation [3]) , where 
Tdust - 45 - 75 K for OHMs and ~ 40 - 120 K for non- 
masing galaxies. 

In contrast, the best fits for both the dust shell thick- 
ness (Y) and optical depth (tv) are markedly different for 
OHMs and non-masing galaxies. The mean Y for non- 
masing galaxies (770 ±320) is nearly twice as thick as the 
mean value for OHMs (350 ± 260), although within the 
large scatter on both parameters. Rather than a physi- 
cal difference in the shell thicknesses, however, we reit- 
erate that this is a likely consequence of the non-masing 
galaxies being better fit by clumpy models (Figure [7J 
and thus a fundamentally different geometry. The dust 
optical depths for the OHMs have a broad distribution 
of Ty between and 450, with more than 50% having 
tv > 350 and a mean of 300. With the exception of 
a single galaxy with t v = 380 (IRAS 06538+4628), all 
non-masing galaxies have Ty < 100 and a median of less 
than 50. 

Overall, fits to the IRS from radiative transfer models 
show a marked difference between the dust geometries 
of masing and non-masing galaxies. Non-masing galax- 
ies are better fit by models with lower optical depth and 
slightly cooler dust temperatures - based on Figure [7J 
this may be due to dust in the form of a clumpy, obscur- 
ing torus (and possibly a partially visible AGN). OHMs 
are almost all well fit by a smooth screen of dust with a 
thinner shell, but with much higher optical depths. We 
emphasize that this "smoothness" is in the context of the 
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Figure 8. Distribution of the dust geometry parameters from DUSTY for the best fits to the IRS spectra. We modeled: (top left) the 
thickness of the dust shell (Y = Router /Rinner), the radial power-law index of the dust density q (top right), the total optical depth Ty 
(bottom left), and the dust temperature at the outer edge (bottom right). OHMs are in red and non-masing galaxies in blue, with the mean 
values for each samples indicated with dashed lines. 



entire nuclear region (and possibly beyond) of the merg- 
ing galaxies. Smaller overdensities within that smooth 
framework are likely sites of star formation and are nec- 
essary to provide the cloud-cloud overlap that produces 
an OHM. 

5. STATISTICAL COMPARISONS 

The primary goal in observing non-masing ULIRGs 
was to directly compare the samples and identify dif- 
ferences that could be triggers of the OHM. Here, we 
present statistical tests comparing the mid-IR and radio 
properties of both samples. 

5.1. Rank correlations of IR and radio properties 

A first-order method of comparing the samples is cor- 
relation between physical parameters. We computed the 
Spearman's rank correlation coefficient (p) for a range 
of properties from Paper I both in the mid-IR and ra- 
dio regimes. For the mid-IR, we tested relationships be- 
tween PAH 6.2 and 11.3 /mi EW and luminosities, con- 
tinuum spectral indices, silicate absorption depths, grey- 
body Tj and constraints on the geometry from the 
DUSTY best-fit models (Table [J). The radio properties 
we explored included the continuum power at 1420 MHz 
(-Pl42o), integrated Lqh (Paper I) and the peak OHM 



flux density at 16 67 MHz (Pi667! iDarling fc Giovanellil 
120001 12001L I2002af ). Results for the Spearman's p tests 
are shown in Table [3] We omit several parameters 
from this table that were measured, but showed no 
significant correlations; these included the mid-IR fine- 
structure line ratios ([Ne III]/[Ne II], [O IV] /[Ne II], and 
[Ne V]/[Ne II]), H2 temperature and gas mass, depth of 
the 6 /xm water ice feature, and the OH hyperfine ratio 

Rh = F 1667/^1665- 

Several of the correlations in Tabic [3] with high signif- 
icance reflect well-known physical relationships; for ex- 
ample, the correlation between P1420 and Lfir. The 
relation of the spectral index 0:30-20 to ' t is ex- 
pected since 0-30-20 samples the Wien side of a black- 
body peaking near 60 fim. The silicate depths at 9.7 and 
18 jxm are also correlated, as expected from the results 
of the DUSTY models. 

The dust temperature and shell thickness Y from the 
DUSTY best-fit models showed a strong anti-correlation 
in OHMs. This picture fits with a smooth dust screen 
enveloping a central source of illumination - thicker shells 
absorb more energy near the inner boundary, resulting in 
a cooler Tdust near the outer boundary. Y -Tdust also was 
the only correlation coefficient that showed significant 
differences between the masing and non-masing samples 
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Table 3 

Spearman rank correlation z-scores for OHMs 





ai5_s 030—20 
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Note. — The Spearman's rank correlation tests statistical dependence of two variables on being monotonic functions, without assuming 
linearity. The s-scores in this table represent the number of standard deviations by which the correlation differs from the null hypothesis 
of no statistical dependence. Correlations higher than 4cr arc in boldface. 



Table 4 

Kolmogorov-Smirnov tests for OHMs and non-masing galaxies 
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Y, and q are 


the best fit of 



IRS spectra to DUSTY models. Results greater than Act significance 
are in boldface. 



{pOHM = —0.97, pnon = —0.64). We attribute this to 
the evidence that non-masing galaxies are poorly fit by 
DUSTY and likely favor a clumpy geometry (Figure [7]). 

For the OH maser emission itself, the only strong cor- 
relation observed was between the peak OHM power and 
the integrated OHM luminosity. No Spearman correla- 
tions between an OHM parameter (Loh, -P1667) and the 
IRS data were found with > 4cr significance, despite the 
fact that several IR features related to dust geometry 
revealed a clear separation of loci between OHMs and 
non-masing galaxies (see Figures l6l and [T"2|) . While these 
parameters have a clear effect on the existence of the 
OHM, the lack of correlation suggests either that specific 
line properties are not well tracked by these parameters, 
or that no single trigger among these is responsible for 
megamaser production. Alternatively, the OHM might 
be the results of stochastic amplification of small-scale 
conditions, with masing simply becoming more common 
when conditions are favorable. 

5.2. Statistical differences between the samples 



A second method of analyzing statistical differences 
between OHMs and non-masing galaxies is the two-sided 
Kolmogorov-Smirnov (K-S) test, which tests the null hy- 
pothesis that the two samples come from the same par- 
ent distribution. Selected results from the K-S tests are 
given in Table 21 where Dks is the maximum separation 
between the scaled cumulative distribution functions and 
N„ is the number of standard deviations by which Dks 
differs from the null hypothesis (ie, the significance of the 
result). We also give the mean values and la standard 
deviations of the properties for both samples. As with 
Table |31 this omits numerous mid-IR and radio proper- 
ties on which we performed K-S tests, but which showed 
no significant difference. 

The majority of the data show K-S results consistent 
with origins from the same distribution; the exceptions 
all relate to dust properties of the galaxies. Two of these 
are quantities directly measured from IRS data: the 30- 
20 /im spectral index and the 9.7 /im silicate depth. The 
remaining significant parameters describe the dust envi- 
ronment modeled by DUSTY: Y, tv, and Tdust all sup- 
port a fundamentally different distribution of the silicate 
dust for the two samples at the 4c level. 

We extended this analysis by refining the greybody 
Tdust measured with Equation[3J which shows only a mild 
significance in the original K-S test (2a). The continuum 
slopes are largely determined by the amount of dust in 
various temperature regimes; however, (^30-20 shows a 
strong difference while ai5_6 (which samples hotter dust) 
does not. This may indicate that only dust in certain 
temperature regimes (ie, the ~ 50 — 100 K region sampled 
by 20-30 fira continuum) is important in triggering OHM 
emission. We tested this by restricting the fit of our 
dust temperature only to data from 20-30 fira, where 
the Wien approximation applies for typical ULIRG dust 
temperatures. In this case, the K-S test yields a much 
higher and statistically significant (6<r) difference for the 
modified dust temperature (730-20) between the OHMs 
and non-masing galaxies. 

K-S tests were also used to quantify the differences seen 
between the samples on the fits with the DUSTY code. 
In particular, the optical depth jy showed a ha difference 
between the two samples, with the typical OHM having 
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Ty a factor of several above a non-masing galaxy (and 
consistent with Dks for the 9.7 /xm feature). The other 
DUSTY parameters show moderate significance (3-5<r); 
however, the results for different distributions of Y likely 
come from a clumpy geometry, rather than a true in- 
crease in the shell thickness ( §4.5[) . Given that the optical 
depth and dust temperature are not independent param- 
eters in DUSTY, all results from the K-S tests strongly 
indicate that the temperature/optical depth of the dust 
(which depends on its geometry) is a key factor in trig- 
gering an OHM. 

We extended the K-S tests by performing a series of 
survival analyses on the same data. Survival analysis 
is particularly suited for flux-limited samples because it 
properly tre ats upper limits for features not detected i n 
all galaxies (|Feigelson fc Nelsodll985t llsobe et al.lll986l) 
We u sed the ASURV package in IRAF (jLavallev et ahl 
1992), which includes the Gehan's generalized Wilcoxon, 
logrank, Peto & Peto, and Peto & Prentice tests. Run- 
ning survival analysis on all measured mid-IR features 
(including atomic and molecular line emission, hydro- 
carbon and gas-phase absorption, PAH, dust and con- 
tinuum features) gave similar results to the K-S tests; 
no parameter showed significant differences between the 
two samples with the exceptions of «3o_2o and S9.7. All 
tests yielded statistically significant differences for these 
features, with an mean significance of 6a for S9.7 and 5<r 
for a 3 o_20- 

Importantly, the results of our survival analysis also 
discount the possibility that other mid-IR features are 
directly related to the presence of an OHM. In particular, 
we detected absorption from hydrocarbons (HACs), gas- 
phase molecules (C2H2, HCN, and CO2) and crystalline 
silicates almost exclusively in the OHM sample (Paper I). 
Using the upper limits on the absorption features in our 
survival analyses, however, we cannot confirm that the 
lack of detections in the non-masing sample is significant 
above the 3c level for any of these parameters. This 
is largely due to lack of sensitivity in the non-masing 
galaxies, since the limits on non-detections are of similar 
magnitudes to the detected absorption in many galaxies. 

6. COMPARING OBSERVATIONS AND THEORY 

Importantly, the IRS data can explore the physics 
of OHMs by testing the predictions of maser pumping 
models. The most recent and comple t e pum ping cal- 
culations come from lLockett fe Elitzun (2008, hereafter 
LE08). The model assumes a slab geometry and uses the 
escape probability method to solve for the level popula- 
tions of the OH molecule. Given assumptions on the 
physical conditions in the masing regions, the overall 
strength of the OHM ( if any) can b e pre dicted. The 
clumpy OHM model of iParra et al.1 (|2005l ) shows that 
the maser optical depth depends most strongly on the 
dust temperature and optical depth. Since both these 
parameters can be estimated from IRS data, our sample 
offers the first opportunity for testing such a model on a 
large number of galaxies. 

The strength of the OHM in the LE08 model is param- 
eterized as the optical depth in the OH line (Tig67, which 
becomes more negative for higher maser gain). To com- 
pare this to observations, we use the line-to-continuum 
ratio to estimate the apparent observed OH optical depth 
(Table P: 



app , ( S1420 + 5*1667 \ ( a \ 

T 1667 = " ln c ' ( 4 ) 

\ ^1420 / 

Here 5*1667 is the pe ak flux density of the OHM at 
1667 M Hz (taken from lDarling fc Giovanelll[2000l l200l 
l2002a|) and 5i42o is the flux density of the radio con- 
tinuum at 1420 MH z from the NRAO VLA Sky Survey 
(|Condon et al.lll998l) . 

The largest uncertainty in Equation's that it assumes 
an OH filling factor of 1; however, VLBI maps of OHM 
galaxies show the OH emission to have both diffuse and 
compact components. In addition, the 1420 MHz radio 
continuum comes from a much larger physical area than 
the OH emission. Therefore, t^ 7 will be a weaker limit 
to the true 1667 MHz optical depth. VLBI observations 
have ma pped the OHM em i ssion for a handful of nearby 
galaxies (lYates et alJfeOOOt iPihlstrom et aLlfeOOlL 120051: 
Klockner et al. 1 120031: ILonsdale et all 120031: IRovilos et al l 
20031: IKlockner fc Baanl 120041: IRichards et al.l 120051: 
Momi ian et al.ll2006[ ). showing that the difference in ap- 
parent optical depth between the entire galaxy and the 
brightest individual maser spots varies by as much as 
At ~ 1 — 4. Furthermore, the gain for individual maser 
spots with cloud-cloud overlap can be as high as several 
hund red (eg, III Zw 35; Diamon d et al.1119 99: Parr a et al.l 
2005), compared to the diffuse background. Since high- 
resolution OH maps do not exist for the vast majority of 
the IRS galaxies, however, we use t^q^ while remaining 
mindful of the above caveats. 

We present tests of the LE08 model parameters using 
two different techniques: one method estimating Tdust 
and Ty directly from the IRS spectra, and a second using 
parameters extracted from the models fit to the IRS data 
using the DUSTY code. 

6.1. Testing the LE08 pumping model with Tdust and Ty 
from IRS data 

The pumping flux of the OHM in the LE08 model de- 
pends most strongly on the pumping flux, which is con- 
trolled by two factors: the self-absorption of the dust 
(depending on Ty) and the Planck function (depending 
on Tdust)- The 9.7 [im silicate feature in the IRS data 
can be used t o estimate the total Ty using the Galactic 
calibration of C|Roche fc Aitkenlll984lh 

t v = (17.0 ± 1.4) x T9.7. (5) 

Secondly, the dust temperature in the masing region 
is estimated from the greybody fit to the IRS and IRAS 
data ( ^4.41 Table [J). Once both parameters for a galaxy 
are estimated, we plot the data on the contours of LE08- 
predicted Ti667 emission (Figure The color of the 
OHM symbols shows their t^J^ on the same scale as the 
LE08 contours. 

According to the standard LE08 model, the most lumi- 
nous OHMs are expected to have Tdust ~ 90 — 150 K and 
tv of a few tens. None of the OHMs with t^q^ < -3.0 
were located near this region; in addition, the parame- 
ters for virtually all observed OHMs lie well away from 
the highest predicted T1667 in Figure [SJ with Tdust = 
40 - 100 K and Ty = 10 - 50. Roughly 15% of the con- 
firmed OHMs have predicted Ti667 > —0.5, which would 
predict almost no masing activity and would fall well be- 
low the limit for inclusion in our sample. 
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Figure 9. Apparent optical depth of the OHM emission (filled 
diamonds) as a function of Ty (Equation^ and dust temperature 
(Table [T}, based on greybody photometry and the 9.7 fim silicate 
depth. Crosses show temperatures and depths of the non-masing 
sample. The error bars in the top left are the average systematic 
uncertainties (or = 20 K, a T = 5). Contours are a zoomed-in 
region of the OHM model of Lockctt k. Elitzur ( 2008). 

To assess the overall fit of the LE08 model, the er- 
rors on our estima t es of Ty and Tdust must be quanti- 
fied. lYun fc Carillil ( 2002( 1 found that accurate tempera- 
ture fits using Eqn. [3] required much higher photometric 
sampling than we possess for this sample, including ra- 
dio and sub-mm data. Oth er SED models (|Fraver et al.l 
119991: iDunne fc Ealesll200ll ) typically fit two- and three- 
component models with differences in dust temperatures 
ranging up to 100 K. As a result, we estimate the av- 
erage <JT dust for each galaxy to be ~ 20 K. This is quite 
high, but is mitigated somewhat by the number of galax- 
ies in our sample. The mean error on the extinction is 
estimated as a Tv — 5, based on the silicate measurement 
technique, the possibility of saturation, and calibration 
in Eqn. [5J Average error bars are shown in the upper 
right corner of Fig. |9] 

The large uncertainties result in considerable scatter 
in the predicted OHM strength, with Ati667 as high as 
1-2 depending on the local gradient of the LE08 model. 
Many of the OHMs lie near contours where Ti667 is a sen- 
sitive function of T dust ; a shift of ~ 5 K could result in 
a change of up to Ari667 = 0.5, while at the same time 
being relatively insensitive to Ty. We show the distri- 
bution of the difference between t^J^ and the predicted 
ri667 from the LE08 model in Figure [TU] The measured 
OH strengths are on average weaker than those predicted 
by the model ((Ari 6 67) = — 0.8); this is consistent, how- 
ever, with the lower bound on rfggy from the OH filling 
factor. The Xreduced f° r the model using this data is 29.1, 
which rejects a correlation hypothesis at the 5<r level. 

While the agreement for individual OHMs is not 
strong, our data is consistent with other predictions of 
the LE08 model. Based on pumping calculations, they 
show that ULIRGs must have a dust temperature greater 
than 45 K in order to achieve population inversion; cooler 
temperatures move the peak of the blackbody too far 
from the main pumping lines to support the necessary 
pumping flux. 90% of the OHMs have T dust > 45 K, 
with the coolest OHM at 37 K; uncertainties of ~ 20 K 
mean that the dust temperatures for all OHMs are fully 
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Figure 10. Distributio n of the difference in Figu re [9] between the 
predicted ri667 from the Lockctt & Elitzur (2008) model and the 
apparent rj'gg'y measured from radio data. The dotted line shows 
the mean of the distribution at Atis67 = —0.8, showing that the 
LE08 model tends to overpredict the strength of the maser. 

consistent with this predicted lower limit. 

Figure [HI also displays the dust parameters for the non- 
masing galaxies. Half of the non-masing galaxies have 
predicted OH luminosities consistent with little masing 
to none at all (such that |ti667| is small), and the Tig67 
predicted by the LE08 model would lie below our de- 
tection threshold of Loh < 10 23 L Q for almost the en- 
tire sample. Based on the LE08 model and IRS data, 
only a single non-masing galaxy (IRAS 23498+2423, 
in the far lower right of Fig. ^ would have been ex- 
pected to show strong megamaser emission. Interest- 
ingly, IRAS 23498+2423 was the object with the highest 
upper limit on maser emission in the non-masing sam- 
ple, with Loh < 2.25. We re-observed this galaxy at the 
Arecibo Observatorj0 in October 2009 to test the LE08 
prediction; no detection of OH was made, confirming an 
upper limit of Lqh < 10 2 27 L r . 



JQ. 



6.2. 



Testing the LE08 pumping model with fits to IRS 
data from DUSTY 

Our second approach for testing the LE08 model cal- 
culates T^st and Ty from the best fits to the IRS data 
with DUSTY (gj3 Table E|. As discussed in ^EM these 
fits give dust temperatures similar to those from the grey- 
body fit, but with optical depths much higher than those 
calculated using only the 9.7 /^m feature, which can be 
saturated in ULIRGs. The Tdust used here is the value 
at the shell's outer edge; since the radial temperature 
profile of the dust is very steep close to the center and 
shallow at the edges (changing by only a few tens of K 
over the outer half of the shell), this temperature repre- 
sents the bulk of the dust mass and is likely a reasonable 
approximation for conditions in the masing regions. 

Figure QT] shows the LE08 predictions with data from 
the DUSTY best fits to the OHMs and non-masing galax- 
ies. We note that LE08 contours are not complete at 
Ty > 300, and that the apparent horizontal feature at 

6 The Arecibo Observatory is part of the National Astronomy 
and Ionosphere Center, which is operated by Cornell University 
under a cooperative agreement with the National Science Founda- 
tion. 
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Figure 11. Apparent optical depth of the OHM emission (filled 
diamonds) as a function of Ty and dust temperature based on 
fits from DUSTY. Black crosses are the non-masing galaxies; the 
error bar in the upper right corner shows the average uncertainty 
(o> ~ 20, ~ 20 K). Contours show the predicted maser strength 
for the LE08 model, which are not complete at ry > 300. Due 
to the coarse gridding of DUSTY models, several galaxies have 
overlapping points on this plot (eg, five galaxies have best fits of 
T dust = 70 K and = 410). 

Ty = 410 is likely an artifact of gridding in the code. 
The distribution of the galaxies is very different from 
that in Figure EJ the OHMs occupy a much larger range 
in optical depth, increasing Ty by an order of magnitude. 
Two distinct loci are visible; the lower left corner con- 
tains the majority of the non-masing galaxies and several 
OHM with T dust ~ 40 - 60 K and t v < 100. The LE08 
model predicts that these galaxies would show little to 
no maser emission. The second group is almost exclu- 
sively composed of OHMs, with warmer temperatures 
(T dust ~ 60 - 80 K) and t v of several hundred. LE08 
predicts a range of Ti667 for these galaxies from —2.0 to 
0. Two non-masing galaxies with warm temperatures do 
not seem associated with either group. 

The overlapping region of OHMs and non-masing 
galaxies with cooler dust temperatures and Ty < 150 
is of particular interest. The dust parameters for these 
galaxies lie well away from the predicted tiqq? peak and 
are close to the minimum predicted inversion tempera- 
ture of 45 K. The overlapping populations at this locus 
include at least one powerful gigamaser, and imply that 
there must be some triggering factor for an OHM be- 
yond Td us t and Tv- This is supported by the fact that all 
OHMs in this region also lie on the horizontal branch of 
the fork diagram (Fig. [6]), classifying them as likely AGN 
hosts. The T dust from DUSTY for all of these OHMs 
is also significantly cooler than that measured with the 
greybody method; a buried AGN might thus be better 
fit with a multi-temperature model. The connection be- 
tween an AGN and OHM suppression, however, is not 
clear; it could represent a different dust geometry not 
conducive to cloud-cloud overlap, or signal a more ad- 
vanced stage in the galaxy merger, thus putting a limit 
on the effective lifetime (and thus observability) of the 
OHM. 

Since the contours in Figure [TT] are incomplete, we can- 
not fully measure the goodness-of-fit in a method similar 
to Figure [TU] Qualitatively, the model makes good pre- 



dictions for the dust parameters for almost all of the non- 
masing galaxies. The exceptions are IRAS 11119+3257 
(log rf 6 r 6 e 7 d ~ -3.0) and IRAS 06538+4628 (log rf 6 r 6 e 7 d ~ 
— 1). The former is the only object whose predicted emis- 
sion lies well above its observational limits on Loh- This 
galaxy is known to show an exc eptional radio exces s 
as measured by its q-parameter (jCondon et al.l 119911 ) , 
with its value of q = 1.2 3 falling well below the mea n 
value for OHMs found bv lDarling k Giovanelhl (|2002a[ >. 
Such an excess commonly indicates that the galaxy hosts 
an AGN; this is supported by the low 6.2 PAH EW 
of IRAS 11119+3257 and its position on the fork dia- 
gram (Fig. El). IRAS 06538+4628 is the only non-masing 
galaxy in our sample for which dust temperatures might 
be too warm to support a strong population inversion. 

There are several reasons why the Tdust and Ty fit from 
DUSTY might differ from methods used in i )6.1l Ty for 
deeply embedded galaxies can suffer from saturation in 
the 9.7 /im feature. Using data from the full dust profile 
(as we do in DUSTY) samples the broader wings of the 
feature it advances up the curve of growth. DUSTY also 
samples the SED at a much higher resolution than the 
greybody fit, albeit in a more limited wavelength regime. 
Finally, the conversion from £9.7 to Ty is based on a 
Galactic calibration; it is not known how dust composi- 
tion might be different in the Milky Way and ULIRGs, 
for example. On the other hand, the DUSTY models re- 
quire the assumption of a specific geometry which may 
not be appropriate (see Fig. [7]) for non-masing galaxies 
and the dustiest OHMs. While neither method is without 
drawbacks, we believe both to have at least some phys- 
ical merit (and are encouraged by the fact that Tdust is 
mostly consistent). 

6.3. Predictions and future observations 

Overall, comparing the IRS data to predictions from 
the LE08 model yielded mixed results. Using parameters 
from the DUSTY code, the LE08 model correctly pre- 
dicted that most non-masing galaxies should have cool 
dust temperatures and low optical depths; however, sev- 
eral megamasers also have Tdust that would be too cool 
for inversion under this model. The observed T"gg* 7 for in- 
dividual sources shows a great deal of scatter; however, 
this is dominated by observational uncertainties in the 
OH filling factor. Both estimates are consistent with the 
LE08 claim that a minimum dust temperature of 45 K is 
required for maser action; within uncertainties, all OHMs 
have Tdust above this value. Based on results from the 
feature-feature diagram, we suggest that future pump- 
ing models include both clumpy and smooth shell dust 
geometries; treatment of OH kinematics might also be 
necessary to model individual sources in more detail. 

Interestingly, the OHM luminosity (which typically 
depends strongly on the total linewidth) does not ap- 
pear to be a strong function of the currently observable 
global host properties. An improved test of the L<jh- 
host galaxy relationship could use VLBI maps of OHM 
galaxies to constrain the true gain in individual clouds to 
determine the filling factor, and then compare these re- 
sults to spatially resolved IR data in the same regions to 
measure Tdust and Ty . If the parameters for OHM pro- 
duction can be fine-tuned based on size scales of ~ 100 pc 
for nearby galaxies, this will greatly assist in comparisons 
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Figure 12. Peak silicate depth at 9.7 /im vs. the spectral slope 
between 20 and 30 fim for OHMs (circles) and non-masing galaxies 
(crosses). For OHMs, the symbol size is proportional to log Lqjj. 
The dashed line shows the rough separation between the loci of 
OHMs and non-masing galaxies. Preselecting OHM hosts based 
on a «3o— 20 — S9.7 cut may be a powerful technique for future 
megamaser surveys. 

of galaxy-wide SEDs for OHMs to non-masing ULIRGs 
at much greater distances. 

The LE08 model depends on a number of other ISM 
properties, some of which can be further constrained by 
the IRS data. These include the ortho-para ratio of H2, 
which affects collision rates and thermalization of the 
gas. LE08 assumes a constant ortho-para ratio of 3; the 
IRS data show that this is only valid for 4/9 ULIRGs 
for which the ratio can be constrained (and can be as 
low as 0.5). The OH column densities measured using 
the 34.6 /im transition lie in the range Noh = 1 — 3 x 
10 17 cm -2 ; this is roughly a factor of two higher than 
the standard value assumed in the LE08 model. 

Results from our IRS data can also narrow poten- 
tial searches for new OHMs, especially at higher red- 
shifts. The most d istant OHM known lies at z = 0.265 
(jBaan et al.l fl9 92) . Since OHMs are associated with 
merging galaxies, which are most plentiful between z ~ 
1 — 3, we expect a higher spatial density of OHMs in 
the early universe (|Darling fc Giovanellil l2002bh . While 
surveys for more distant OHMs are restricted both by 
sensitivity constraints and low-frequency RFI, a signifi- 
cant obstacle has been identification of a suitable target 
sample of host galaxies. Based on the IRS data, we sug- 
gest that future OHM surveys target galaxies with dust 
peaks near X res t = 53 /im, steep 30 — 20 ^m slopes, deep 
dust absorption, and that do not show evidence of host- 
ing an AGN. Figure [12] shows how the combination of 
&30-20 and S9.7 can clearly separate almost all OHMs 
from non-masing galaxies in the mid-infrared. This may 
be a valuable tool in future searches for OHMs; for galax- 
ies in which low-resolution IR spectroscopy is available, 
pre-selecting OHM candidates based on these diagnos- 
tics should have a success rate far in excess of blindly 
selecting ULIRGs from the field. The growing number of 
sub-millimeter galaxy catalogs and multiwavelength deep 
fields offer excellent opportunities in the near future for 
such surveys. 

Finally, new observatories are offering opportunities 
for completing the OHM picture. There are only a few 
OHMs in which the important 53 /im transition has been 



measured (eg, iHe fc Ch en 2004); this could be poten- 
tially obserbed in larger numbers of galaxies with SOFIA. 
The Herschel observatory can also supplement the IRS 
by measuring the 79 and 119 um OH transit ions in large 
samples of ULIRGs (eg, iFischer et alll2010[ ). If models 
such as LE08 can be refined based on direct OH mea- 
surements, then the pumping efficiencies for the mega- 
maser could be evaluated for a statistically significant 
sample. Photometric measurements from these instru- 
ments and from JWST will also generate SEDs with 
much broader spectral coverage, increasing our knowl- 
edge of the physics needed for radiative transfer models 
and modeling the OHM environment. 

7. CONCLUSIONS 

We present results from the Spitzer Infrared Spectro- 
graph comparing the mid-IR properties of OH mega- 
maser hosts to galaxies with confirmed upper limits on 
the megamaser emission. No significant differences be- 
tween the samples were found for the average excitation 
states, line velocities, or star formation rates. 10-25% 
of the OHMs show clear evidence in the mid-IR for an 
AGN, significantly lower than previous optical and radio 
studies which placed the AGN fraction of OHMs between 
30-70%. In non-masing ULIRGs, between 50-60% of the 
galaxies have mid-IR evidence for an AGN. 

Fits of radiative transfer models to the IRS spectra 
with the DUSTY code show that OHMs have warmer 
Tdust and deeper silicate absorption associated with a 
smooth, thick dust shell surrounding the nucleus. This 
implies the presence of a large dust reservoir in OHMs 
with a smooth geometry and temperatures from ~ 50 — 
100 K. Non-masing galaxies show weaker dust absorp- 
tion, shallower mid-IR continuum, and cooler dust (by 
~ 10 K) than the typical OHM host. The relative 
strength of silicate features in non-masing galaxies sug- 
gests that they are best fit with a clumpy dust geometry. 

We used IRS data to evaluate predictions from the 
OH pumping model of lLockett fc Elitzurl (|2008l) . the first 
direct test of OHM production using observed proper- 
ties of the host galaxies. The dust opacities for OHMs 
derived from the best-fit DUSTY models suggest that 
much higher opacities (ry ~ 100 — 400) are necessary 
for OHM production. All the IRS data are consistent 
with the LE08 claim that a minimum Tdust — 45 K is 
required for maser action. Limits on the OH emission 
for most non-masing galaxies are predicted by the LE08 
model, based on their comparatively cool dust tempera- 
tures (Tdust < 60 K) and low dust opacity (ry < 100). 
Finally, the IRS data constrain several parameters neces- 
sary to develop future OHM pumping models, including 
the dust optical depth, temperature, and overall geome- 
try. 

For the first time, we present spectral diagnostics that 
can distinguish OHMs from non-masing galaxies based 
on their host galaxy properties, the clearest of which is 
the Sg.7 — 030-20 relation (Figure[T2]). These parameters 
can be relatively easily measured in low-resolution spec- 
tra, and may signify a powerful method for preselecting 
OHM candidates for follow-up radio surveys at higher 
redshift. 

This work is based on observations made with the 
Spitzer Space Telescope, which is operated by the Jet 
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Propulsion Laboratory (JPL) and the California Insti- 
tute of Technology, under NASA contract 1407. Addi- 
tional observations were also taken at the Arecibo Obser- 
vatory (NAIC/Cornell). We made extensive use of the 
NASA/IPAC Extragalactic Database (NED) which is op- 
erated by JPL and Caltech under contract with NASA. 
Many thanks to are due to J. Stocke for jump-starting 
work with Spitzer, D. Dale for discussions about mid-IR 
line ratios, P. Lockett and M. Elitzur for sharing their 
radiative transfer models, and the Spitzer Science Cen- 
ter for hosting KW and JD while collaborating on data 
analysis. VC acknowledges partial support from the EU 
ToK grant 39965 and FP7-REGPOT 206469. 
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